SLAMseq resolves the kinetics of maternal and zygotic gene expression during early zebrafish embryogenesis.

The maternal-to-zygotic transition (MZT) is a key developmental process in metazoan embryos that involves the activation of zygotic transcription (ZGA) and degradation of maternal transcripts. We employed metabolic mRNA sequencing (SLAMseq) to deconvolute the compound embryonic transcriptome in zebrafish. While mitochondrial zygotic transcripts prevail prior to MZT, we uncover the spurious transcription of hundreds of short and intron-poor genes as early as the 2-cell stage. Upon ZGA, most zygotic transcripts originate from thousands of maternal-zygotic (MZ) genes that are transcribed at rates comparable to those of hundreds of purely zygotic genes and replenish maternal mRNAs at distinct timescales. Rapid replacement of MZ transcripts involves transcript decay features unrelated to major maternal degradation pathways and promotes de novo synthesis of the core gene expression machinery by increasing poly(A)-tail length and translation efficiency. SLAMseq hence provides insights into the timescales, molecular features, and regulation of MZT during zebrafish embryogenesis.

Hijacking of transcriptional condensates by endogenous retroviruses.

Most endogenous retroviruses (ERVs) in mammals are incapable of retrotransposition; therefore, why ERV derepression is associated with lethality during early development has been a mystery. Here, we report that rapid and selective degradation of the heterochromatin adapter protein TRIM28 triggers dissociation of transcriptional condensates from loci encoding super-enhancer (SE)-driven pluripotency genes and their association with transcribed ERV loci in murine embryonic stem cells. Knockdown of ERV RNAs or forced expression of SE-enriched transcription factors rescued condensate localization at SEs in TRIM28-degraded cells. In a biochemical reconstitution system, ERV RNA facilitated partitioning of RNA polymerase II and the Mediator coactivator into phase-separated droplets. In TRIM28 knockout mouse embryos, single-cell RNA-seq analysis revealed specific depletion of pluripotent lineages. We propose that coding and noncoding nascent RNAs, including those produced by retrotransposons, may facilitate 'hijacking' of transcriptional condensates in various developmental and disease contexts.

NMD is required for timely cell fate transitions by fine-tuning gene expression and regulating translation.

Cell fate transitions depend on balanced rewiring of transcription and translation programs to mediate ordered developmental progression. Components of the nonsense-mediated mRNA decay (NMD) pathway have been implicated in regulating embryonic stem cell (ESC) differentiation, but the exact mechanism is unclear. Here we show that NMD controls expression levels of the translation initiation factor Eif4a2 and its premature termination codon-encoding isoform (Eif4a2PTC ). NMD deficiency leads to translation of the truncated eIF4A2PTC protein. eIF4A2PTC elicits increased mTORC1 activity and translation rates and causes differentiation delays. This establishes a previously unknown feedback loop between NMD and translation initiation. Furthermore, our results show a clear hierarchy in the severity of target deregulation and differentiation phenotypes between NMD effector KOs (Smg5 KO > Smg6 KO > Smg7 KO), which highlights heterodimer-independent functions for SMG5 and SMG7. Together, our findings expose an intricate link between mRNA homeostasis and mTORC1 activity that must be maintained for normal dynamics of cell state transitions.

Transcriptome-Wide Profiling of RNA Stability.

Gene expression is controlled at multiple levels, including RNA transcription and turnover. But determining the relative contributions of RNA biogenesis and decay to the steady-state abundance of cellular transcripts remains challenging because conventional transcriptomics approaches do not provide the temporal resolution to derive the kinetic parameters underlying steady-state gene expression.Here, we describe a protocol that combines metabolic RNA labeling by 4-thiouridine with chemical nucleoside conversion and whole-transcriptome sequencing followed by bioinformatics analysis to determine RNA stability in cultured cells at a genomic scale. Time-resolved transcriptomics by thiol (SH)-linked alkylation for the metabolic sequencing of RNA (SLAMseq) provides accurate information on transcript half-lives across annotated features in the genome, including by-products of transcription, such as introns. We provide a step-by-step instruction for time-resolved transcriptomics, which enhances traditional RNA sequencing protocols to acquire the temporal resolution required to directly measure the cellular kinetics of RNA turnover under physiological conditions.

Systematic refinement of gene annotations by parsing mRNA 3' end sequencing datasets.

Alternative cleavage and polyadenylation generates mRNA 3' isoforms in a cell type-specific manner. Due to finite available RNA sequencing data of organisms with vast cell type complexity, currently available gene annotation resources are incomplete, which poses significant challenges to the comprehensive interpretation and quantification of transcriptomes. In this chapter, we introduce 3'GAmES, a stand-alone computational pipeline for the identification and quantification of novel mRNA 3'end isoforms from 3'mRNA sequencing data. 3'GAmES expands available repositories and improves comprehensive gene-tag counting by cost-effective 3' mRNA sequencing, faithfully mirroring whole-transcriptome RNAseq measurements. By employing R and bash shell scripts (assembled in a Singularity container) 3'GAmES systematically augments cell type-specific 3' ends of RNA polymerase II transcripts and increases the sensitivity of quantitative gene expression profiling by 3' mRNA sequencing. Public access: https://github.com/AmeresLab/3-GAmES.git.

Structure-function analysis of microRNA 3'-end trimming by Nibbler.

Nibbler (Nbr) is a 3'-to-5' exoribonuclease whose catalytic 3'-end trimming activity impacts microRNA (miRNA) and PIWI-interacting RNA (piRNA) biogenesis. Here, we report on structural and functional studies to decipher the contributions of Nbr's N-terminal domain (NTD) and exonucleolytic domain (EXO) in miRNA 3'-end trimming. We have solved the crystal structures of the NTD core and EXO domains of Nbr, both in the apo-state. The NTD-core domain of Aedes aegypti Nbr adopts a HEAT-like repeat scaffold with basic patches constituting an RNA-binding surface exhibiting a preference for binding double-strand RNA (dsRNA) over single-strand RNA (ssRNA). Structure-guided functional assays in Drosophila S2 cells confirmed a principal role of the NTD in exonucleolytic miRNA trimming, which depends on basic surface patches. Gain-of-function experiments revealed a potential role of the NTD in recruiting Nbr to Argonaute-bound small RNA substrates. The EXO domain of A. aegypti and Drosophila melanogaster Nbr adopt a mixed α/ß-scaffold with a deep pocket lined by a DEDDy catalytic cleavage motif. We demonstrate that Nbr's EXO domain exhibits Mn2+-dependent ssRNA-specific 3'-to-5' exoribonuclease activity. Modeling of a 3' terminal Uridine into the catalytic pocket of Nbr EXO indicates that 2'-O-methylation of the 3'-U would result in a steric clash with a tryptophan side chain, suggesting that 2'-O-methylation protects small RNAs from Nbr-mediated trimming. Overall, our data establish that Nbr requires its NTD as a substrate recruitment platform to execute exonucleolytic miRNA maturation, catalyzed by the ribonuclease EXO domain.

Conformation of sister chromatids in the replicated human genome.

The three-dimensional organization of the genome supports regulated gene expression, recombination, DNA repair, and chromosome segregation during mitosis. Chromosome conformation capture (Hi-C)1,2 analysis has revealed a complex genomic landscape of internal chromosomal structures in vertebrate cells3-7, but the identical sequence of sister chromatids has made it difficult to determine how they topologically interact in replicated chromosomes. Here we describe sister-chromatid-sensitive Hi-C (scsHi-C), which is based on labelling of nascent DNA with 4-thio-thymidine and nucleoside conversion chemistry. Genome-wide conformation maps of human chromosomes reveal that sister-chromatid pairs interact most frequently at the boundaries of topologically associating domains (TADs). Continuous loading of a dynamic cohesin pool separates sister-chromatid pairs inside TADs and is required to focus sister-chromatid contacts at TAD boundaries. We identified a subset of TADs that are overall highly paired and are characterized by facultative heterochromatin and insulated topological domains that form separately within individual sister chromatids. The rich pattern of sister-chromatid topologies and our scsHi-C technology will make it possible to investigate how physical interactions between identical DNA molecules contribute to DNA repair, gene expression, chromosome segregation, and potentially other biological processes.

Determining mRNA Stability by Metabolic RNA Labeling and Chemical Nucleoside Conversion.

The varying rates at which mRNAs decay are tightly coordinated with transcriptional changes to shape gene expression during development and disease. But currently available RNA sequencing approaches lack the temporal information to determine the relative contribution of RNA biogenesis, processing and turnover to the establishment of steady-state gene expression profiles.Here, we describe a protocol that combines metabolic RNA labeling with chemical nucleoside conversion by thiol-linked alkylation of 4-thiouridine to determine RNA stability in cultured cells (SLAMseq). When coupled to cost-effective mRNA 3' end sequencing approaches, SLAMseq determines the half-life of polyadenylated transcripts in a global and transcript-specific manner using untargeted or targeted cDNA library preparation protocols.We provide a step-by-step instruction for time-resolved mRNA 3' end sequencing, which augments traditional RNA-seq approaches to acquire the temporal resolution necessary to study the molecular principles that control gene expression.

Publisher Correction: Sequencing cell-type-specific transcriptomes with SLAM-ITseq.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Time-Resolved Small RNA Sequencing Unravels the Molecular Principles of MicroRNA Homeostasis.

Argonaute-bound microRNAs silence mRNA expression in a dynamic and regulated manner to control organismal development, physiology, and disease. We employed metabolic small RNA sequencing for a comprehensive view on intracellular microRNA kinetics in Drosophila. Based on absolute rate of biogenesis and decay, microRNAs rank among the fastest produced and longest-lived cellular transcripts, disposing up to 105 copies per cell at steady-state. Mature microRNAs are produced within minutes, revealing tight intracellular coupling of biogenesis that is selectively disrupted by pre-miRNA-uridylation. Control over Argonaute protein homeostasis generates a kinetic bottleneck that cooperates with non-coding RNA surveillance to ensure faithful microRNA loading. Finally, regulated small RNA decay enables the selective rapid turnover of Ago1-bound microRNAs, but not of Ago2-bound small interfering RNAs (siRNAs), reflecting key differences in the robustness of small RNA silencing pathways. Time-resolved small RNA sequencing opens new experimental avenues to deconvolute the timescales, molecular features, and regulation of small RNA silencing pathways in living cells.

Sequencing cell-type-specific transcriptomes with SLAM-ITseq.

Analysis of cell-type-specific transcriptomes is vital for understanding the biology of tissues and organs in the context of multicellular organisms. In this Protocol Extension, we combine a previously developed cell-type-specific metabolic RNA labeling method (thiouracil (TU) tagging) and a pipeline to detect the labeled transcripts by a novel RNA sequencing (RNA-seq) method, SLAMseq (thiol (SH)-linked alkylation for the metabolic sequencing of RNA). By injecting a uracil analog, 4-thiouracil, into transgenic mice that express cell-type-specific uracil phosphoribosyltransferase (UPRT), an enzyme required for 4-thiouracil incorporation into newly synthesized RNA, only cells expressing UPRT synthesize thiol-containing RNA. Total RNA isolated from a tissue of interest is then sequenced with SLAMseq, which introduces thymine to cytosine (T>C) conversions at the sites of the incorporated 4-thiouracil. The resulting sequencing reads are then mapped with the T>C-aware alignment software, SLAM-DUNK, which allows mapping of reads containing T>C mismatches. The number of T>C conversions per transcript is further analyzed to identify which transcripts are synthesized in the UPRT-expressing cells. Thus, our method, SLAM-ITseq (SLAMseq in tissue), enables cell-specific transcriptomics without laborious FACS-based cell sorting or biochemical isolation of the labeled transcripts used in TU tagging. In the murine tissues we assessed previously, this method identified ~5,000 genes that are expressed in a cell type of interest from the total RNA pool from the tissue. Any laboratory with access to a high-throughput sequencer and high-power computing can adapt this protocol with ease, and the entire pipeline can be completed in <5 d.

Quantification of experimentally induced nucleotide conversions in high-throughput sequencing datasets.

BACKGROUND: Methods to read out naturally occurring or experimentally introduced nucleic acid modifications are emerging as powerful tools to study dynamic cellular processes. The recovery, quantification and interpretation of such events in high-throughput sequencing datasets demands specialized bioinformatics approaches. RESULTS: Here, we present Digital Unmasking of Nucleotide conversions in K-mers (DUNK), a data analysis pipeline enabling the quantification of nucleotide conversions in high-throughput sequencing datasets. We demonstrate using experimentally generated and simulated datasets that DUNK allows constant mapping rates irrespective of nucleotide-conversion rates, promotes the recovery of multimapping reads and employs Single Nucleotide Polymorphism (SNP) masking to uncouple true SNPs from nucleotide conversions to facilitate a robust and sensitive quantification of nucleotide-conversions. As a first application, we implement this strategy as SLAM-DUNK for the analysis of SLAMseq profiles, in which 4-thiouridine-labeled transcripts are detected based on T > C conversions. SLAM-DUNK provides both raw counts of nucleotide-conversion containing reads as well as a base-content and read coverage normalized approach for estimating the fractions of labeled transcripts as readout. CONCLUSION: Beyond providing a readily accessible tool for analyzing SLAMseq and related time-resolved RNA sequencing methods (TimeLapse-seq, TUC-seq), DUNK establishes a broadly applicable strategy for quantifying nucleotide conversions.

Structural basis for acceptor RNA substrate selectivity of the 3' terminal uridylyl transferase Tailor.

Non-templated 3'-uridylation of RNAs has emerged as an important mechanism for regulating the processing, stability and biological function of eukaryotic transcripts. In Drosophila, oligouridine tailing by the terminal uridylyl transferase (TUTase) Tailor of numerous RNAs induces their degradation by the exonuclease Dis3L2, which serves functional roles in RNA surveillance and mirtron RNA biogenesis. Tailor preferentially uridylates RNAs terminating in guanosine or uridine nucleotides but the structural basis underpinning its RNA substrate selectivity is unknown. Here, we report crystal structures of Tailor bound to a donor substrate analog or mono- and oligouridylated RNA products. These structures reveal specific amino acid residues involved in donor and acceptor substrate recognition, and complementary biochemical assays confirm the critical role of an active site arginine in conferring selectivity toward 3'-guanosine terminated RNAs. Notably, conservation of these active site features suggests that other eukaryotic TUTases, including mammalian TUT4 and TUT7, might exhibit similar, hitherto unknown, substrate selectivity. Together, these studies provide critical insights into the specificity of 3'-uridylation in eukaryotic post-transcriptional gene regulation.

Small RNAs Are Trafficked from the Epididymis to Developing Mammalian Sperm.

The biogenesis of the RNA payload of mature sperm is of great interest, because RNAs delivered to the zygote at fertilization can affect early development. Here, we tested the hypothesis that small RNAs are trafficked to mammalian sperm during the process of post-testicular maturation in the epididymis. By characterizing small RNA dynamics during germ cell maturation in mice, we confirm and extend prior observations that sperm undergo a dramatic switch in the RNA payload from piRNAs to tRNA fragments (tRFs) upon exiting the testis and entering the epididymis. Small RNA delivery to sperm could be recapitulated in vitro by incubating testicular spermatozoa with caput epididymosomes. Finally, tissue-specific metabolic labeling of RNAs in intact mice definitively shows that mature sperm carry RNAs that were originally synthesized in the epididymal epithelium. These data demonstrate that soma-germline RNA transfer occurs in male mammals, most likely via vesicular transport from the epididymis to maturing sperm.

Analysis of 3' End Modifications in microRNAs by High-Throughput Sequencing.

MicroRNAs are ~22 nt small, non-coding RNAs that direct posttranscriptional silencing of gene expression to regulate animal development, physiology, and disease. An emerging mechanism that controls the biogenesis of microRNAs is the addition of non-templated nucleotides, predominantly uridine, to the 3' end of precursor-microRNAs, in a process that is commonly referred to as tailing. Here, we describe methods that enable the systematic characterization of tailing events in mature microRNAs and their precursors. We report protocols for untargeted and targeted cDNA library preparation procedures, as exemplified in the context of the model organism Drosophila melanogaster and focusing on precursor-microRNAs. We also refer to a dedicated computational framework for the subsequent analysis of untemplated nucleotide additions in cDNA libraries. The described methods for the systematic characterization of posttranscriptional modifications in gene regulatory small RNAs and their precursors will be instrumental in clarifying regulatory concepts that control posttranscriptional gene silencing.

SLAM-ITseq: sequencing cell type-specific transcriptomes without cell sorting.

Cell type-specific transcriptome analysis is an essential tool for understanding biological processes in which diverse types of cells are involved. Although cell isolation methods such as fluorescence-activated cell sorting (FACS) in combination with transcriptome analysis have widely been used so far, their time-consuming and harsh procedures limit their applications. Here, we report a novel in vivo metabolic RNA sequencing method, SLAM-ITseq, which metabolically labels RNA with 4-thiouracil in a specific cell type in vivo followed by detection through an RNA-seq-based method that specifically distinguishes the thiolated uridine by base conversion. This method has successfully identified the cell type-specific transcriptome in three different tissues: endothelial cells in brain, epithelial cells in intestine and adipocytes in white adipose tissue. As this method does not require isolation of cells or RNA prior to the transcriptomic analysis, SLAM-ITseq provides an easy yet accurate snapshot of the transcriptional state in vivo.

SLAM-seq defines direct gene-regulatory functions of the BRD4-MYC axis.

Defining direct targets of transcription factors and regulatory pathways is key to understanding their roles in physiology and disease. We combined SLAM-seq [thiol(SH)-linked alkylation for the metabolic sequencing of RNA], a method for direct quantification of newly synthesized messenger RNAs (mRNAs), with pharmacological and chemical-genetic perturbation in order to define regulatory functions of two transcriptional hubs in cancer, BRD4 and MYC, and to interrogate direct responses to BET bromodomain inhibitors (BETis). We found that BRD4 acts as general coactivator of RNA polymerase II-dependent transcription, which is broadly repressed upon high-dose BETi treatment. At doses triggering selective effects in leukemia, BETis deregulate a small set of hypersensitive targets including MYC. In contrast to BRD4, MYC primarily acts as a selective transcriptional activator controlling metabolic processes such as ribosome biogenesis and de novo purine synthesis. Our study establishes a simple and scalable strategy to identify direct transcriptional targets of any gene or pathway.

Cell-type specific sequencing of microRNAs from complex animal tissues.

MicroRNAs (miRNAs) play an essential role in the post-transcriptional regulation of animal development and physiology. However, in vivo studies aimed at linking miRNA function to the biology of distinct cell types within complex tissues remain challenging, partly because in vivo miRNA-profiling methods lack cellular resolution. We report microRNome by methylation-dependent sequencing (mime-seq), an in vivo enzymatic small-RNA-tagging approach that enables high-throughput sequencing of tissue- and cell-type-specific miRNAs in animals. The method combines cell-type-specific 3'-terminal 2'-O-methylation of animal miRNAs by a genetically encoded, plant-specific methyltransferase (HEN1), with chemoselective small-RNA cloning and high-throughput sequencing. We show that mime-seq uncovers the miRNomes of specific cells within Caenorhabditis elegans and Drosophila at unprecedented specificity and sensitivity, enabling miRNA profiling with single-cell resolution in whole animals. Mime-seq overcomes current challenges in cell-type-specific small-RNA profiling and provides novel entry points for understanding the function of miRNAs in spatially restricted physiological settings.

Thiol-linked alkylation of RNA to assess expression dynamics.

Gene expression profiling by high-throughput sequencing reveals qualitative and quantitative changes in RNA species at steady state but obscures the intracellular dynamics of RNA transcription, processing and decay. We developed thiol(SH)-linked alkylation for the metabolic sequencing of RNA (SLAM seq), an orthogonal-chemistry-based RNA sequencing technology that detects 4-thiouridine (s4U) incorporation in RNA species at single-nucleotide resolution. In combination with well-established metabolic RNA labeling protocols and coupled to standard, low-input, high-throughput RNA sequencing methods, SLAM seq enabled rapid access to RNA-polymerase-II-dependent gene expression dynamics in the context of total RNA. We validated the method in mouse embryonic stem cells by showing that the RNA-polymerase-II-dependent transcriptional output scaled with Oct4/Sox2/Nanog-defined enhancer activity, and we provide quantitative and mechanistic evidence for transcript-specific RNA turnover mediated by post-transcriptional gene regulatory pathways initiated by microRNAs and N6-methyladenosine. SLAM seq facilitates the dissection of fundamental mechanisms that control gene expression in an accessible, cost-effective and scalable manner.

Genetic and mechanistic diversity of piRNA 3'-end formation.

Small regulatory RNAs guide Argonaute (Ago) proteins in a sequence-specific manner to their targets and therefore have important roles in eukaryotic gene silencing. Of the three small RNA classes, microRNAs and short interfering RNAs are processed from double-stranded precursors into defined 21- to 23-mers by Dicer, an endoribonuclease with intrinsic ruler function. PIWI-interacting RNAs (piRNAs)-the 22-30-nt-long guides for PIWI-clade Ago proteins that silence transposons in animal gonads-are generated independently of Dicer from single-stranded precursors. piRNA 5' ends are defined either by Zucchini, the Drosophila homologue of mitoPLD-a mitochondria-anchored endonuclease, or by piRNA-guided target cleavage. Formation of piRNA 3' ends is poorly understood. Here we report that two genetically and mechanistically distinct pathways generate piRNA 3' ends in Drosophila. The initiating nucleases are either Zucchini or the PIWI-clade proteins Aubergine (Aub) or Ago3. While Zucchini-mediated cleavages directly define mature piRNA 3' ends, Aub/Ago3-mediated cleavages liberate pre-piRNAs that require extensive resection by the 3'-to-5' exoribonuclease Nibbler (Drosophila homologue of Mut-7). The relative activity of these two pathways dictates the extent to which piRNAs are directed to cytoplasmic or nuclear PIWI-clade proteins and thereby sets the balance between post-transcriptional and transcriptional silencing. Notably, loss of both Zucchini and Nibbler reveals a minimal, Argonaute-driven small RNA biogenesis pathway in which piRNA 5' and 3' ends are directly produced by closely spaced Aub/Ago3-mediated cleavage events. Our data reveal a coherent model for piRNA biogenesis, and should aid the mechanistic dissection of the processes that govern piRNA 3'-end formation.